Dissipation-induced symmetry breaking in a driven optical lattice 
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We analyze the atomic dynamics in an ac driven periodic optical potential which is symmetric 
in both time and space. We experimentally demonstrate that in the presence of dissipation the 
symmetry is broken, and a current of atoms through the optical lattice is generated as a result. 
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Brownian motors 0, IE IE El IE IE 01 are devices which 
rectify the random motion of Brownian particles, gener- 
ating in this way a current. Quite recently a large amount 
of work has been devoted to the study of Brownian mo- 
tors @, H EE EE 13 EE ES EE HE Gj; as they seem to 
have important applications in very different areas. On 
one hand, Brownian motors may be the key for the un- 
derstanding of the working principle of molecular motors, 
tiny biological engines which transform the energy pro- 
duced in chemical reactions into unidirectional motion 
along periodic structures which are macroscopically flat 
. On the other hand, the mechanism of rectification of 
fluctuations identified in the study of Brownian motors 
may lead to new electron pumps, and indeed the study of 
solid state devices which implement Brownian motors is 
at present a very active area of research [l8|, lltl l20i |21| . 

The realization of a Brownian motor is a quite subtle 
task, and requires to overcome the limitations imposed 
by the second principle of thermodynamics and to break 
the symmetries which inhibit directed motion. Indeed, as 
the second principle of thermodynamics does not allow 
the appearance of a current in a system at thermal equi- 
librium, Brownian motors are realized by driving Brown- 
ian particles out of equilibrium, as identified in the early 
proposals for flashing p| and rocking [E 01 ratchets. In 
order then to obtain directed motion in the system out 
of equilibrium, relevant symmetries have to be broken. 

Theoretical work 0, 13 clearly identified the symme- 
tries which in the Hamiltonian limit, i.e. in the ab- 
sence of dissipation, inhibit directed motion. However, 
in the pres ence of dissipation the scenario may change 
[Til Il5l.ll6| and it was theoretically shown that the sym- 
metry properties of the system are modified by the pres- 
ence of dissipation and a current can be generated even 
when the symmetries of the Hamiltonian would prevent 
current generation in the Hamiltonian limit |l4 . Il5| . 

In this work we demonstrate experimentally the phe- 
nomenon of dissipation-induced symmetry breaking us- 
ing cold atoms in an ac driven periodic optical potential 
|22j which is symmetric in both time and space. We show 
that in the presence of dissipation the symmetry is bro- 
ken, and a current of atoms through the optical lattice is 
generated as a result. 

Before describing our experimental observations, it 
is essential to introduce the relevant symmetries which 
control the current generation in the Hamiltonian limit 



EE lldl ITU . We consider a particle in a spatially symmet- 
ric periodic potential, periodically rocked by a zero-mean 
force F(t) of period T = 2tt/uj. In this case there are two 
relevant symmetries which need to be examined to deter- 
mine whether current generation is possible. Following 
the notations of Ref. EH we say that F(t) possesses F s 
symmetry if F(t) is symmetric, after some appropriate 
shift: F(t + t) = F(-t + r). Moreover, if F(t) satisfies 
F(t) = —F(t + T/2), we say that F possesses F s h sym- 
metry. The symmetry F s h of the driving implies that 
the system is invariant under the shift-transformation 
S a : (x,p,t) — > (— x, — p, t + T/2), while the symmetry 
F s leads to invariance under the time-reversal transfor- 
mation Sb ■ (x,p,t) — > (x,—p,—t). The invariance of 
the system under any of these two transformations for- 
bids the appearance of a directed current 0, [^. To 
elaborate further, we consider the specific form for F(t): 



F(t) = F Q [A cos cot + B cos(2wt + 4>)\ 



(1) 



The presence of both harmonics (Fq,A, B ^ 0) breaks 
the shift symmetry F s h, independently of the value of 
the relative phase <f>. On the other hand, whether or 
not the F s symmetry is broken depends on the value of 
the phase <fi: for <fi = mt, with n integer, the symmetry 
F s is preserved, while for <j) ^ rm it is broken. There- 
fore for <f> — nir current generation is forbidden in the 
Hamiltonian limit, while for cf> ^= nn is allowed. Per- 
turbative calculations |14j show that the average current 
of particles is, in leading order, proportional to sin 0, in 
agreement with the symmetry considerations discussed 
above. The limit in which the effect of dissipation on the 
symmetries is negligible was already experimentally ex- 
amined in Ref. [23j where the dependence / ~ sin (j> was 
demonstrated. 

Consider now the case of weak, nonzero dissipation. 
The presence of dissipation breaks the invariance under 
time-reversal transformation Sb even if the driving is F s - 
symmetric. Therefore for a biharmonic force of the form 
of Eq. both time-reversal and shift symmetries are 
violated in the presence of dissipation. A current can be 
generated as a result. An interesting issue is how the 
dependence of the current on the phase <f> is modified by 
the presence of dissipation. Calculations done by solving 
the kinetic Boltzmann equation for an ensemble of inter- 
acting particles EE] showed that in the presence of weak 
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dissipation the average current of particles / still shows 
an approximate sinusoidal dependence on the phase </>, 
but a phase lag 4>o appears as a result of dissipation: 

/ ~ sin(0 - 0q) . (2) 

The phase lag 4>o vanishes in the Hamiltonian limit and 
is an increasing function of the relaxation rate [I5f . This 
dependence is consistant with the previous observation 
that weak dissipation breaks the time-reversal symme- 
try of the system and leads to the generation of current 
also for (j> = nn. In other words, the shift </>o is a signa- 
ture of the phenomenon of dissipation-induced symmetry 
breaking, and in our experiment </> will precisely be the 
quantity examined to detect the phenomenon. 

In our experiment we demonstrate the phenomenon 
of dissipation-induced symmetry breaking by using cold 
atoms in an ac driven near-resonant optical lattice |22| . 
This is the same system used previously to demonstrate 
the rectification of fluctuations [24| and to investigate the 
phenomenon of resonant activation |25| . These investi- 
gations clearly showed that near-resonant optical lattices 
represent an ideal model system to investigate phenom- 
ena of statistical physics. In fact, in near-resonant optical 
lattices the laser fields create both a periodic potential 
for the atoms and introduce dissipation. More precisely, 
the interference of the laser fields creates one periodic 
potential for each ground state of the atoms. The laser 
fields also introduce stochastic transitions (optical pump- 
ing processes) between different ground states. This leads 
to damping, an effect named Sisyphus cooling, and a fluc- 
tuating force. As a result of the fluctuations, the atoms 
undergo a random walk through the periodic potential, 
and indeed normal diffusion was observed for atoms in 
dissipative optical lattices pfiq . 

The experimental set-up is the same as the one used 
in our previous work |24l l25| . Cesium atoms are cooled 
and trapped in a 3D optical lattice created by four laser 
beams arranged in the so-called umbrella-like configura- 
tion. One beam (beam 1) propagates in the z-direction; 
the three other beams (beams 2-4) propagate in the op- 
posite direction, arranged along the edges of a triangular 
pyramid having the z direction as axis. For further detail 
on the lattice beams arrangement, as lattice beams' an- 
gles and polarizations, we refer to Ref. j^. A zero-mean 
oscillating force of the form can be applied by phase 
modulating one of the lattice beams. More precisley [24| 
a phase modulation of beam 1 of the form 

a(t) — a n [A cos ujt + — cos(2ujt — (j))] (3) 

will result, in the accelerated frame in which the opti- 
cal lattice is stationary, in an homogeneous force in the 
z direction of the form of Eq. Q with F = mLU 2 a /2k, 
where m is the atomic mass and k the laser field wavevec- 
tor. 

Before describing the experimental results, it is nec- 
essary to analyze theoretically our system to determine 



whether the description of the dissipation-induced sym- 
metry breaking in terms of a phase lag <p derived for 
an ensemble of particles in the presence of collisions 0] 
applies also to the present case of non-interacting atoms, 
with the dissipation associated to the scattering of pho- 
tons. For the sake of simplicity, and to make the analysis 
more transparent, we consider the simplest atom-light 
configuration in which Sisyphus cooling takes place: a 
J g = 1/2 — > J e = 3/2 atomic transition, and a ID 
light configuration consisting of two linearly polarized 
laser fields, counterpropagating and with orthogonal po- 
larizations - the so called lin_Llin configuration 22] . This 
atom-light configuration results in two optical potentials 
U± for the atoms, one for each ground state |±), in phase 
opposition: U± = Uq[— 2 ± cos2fcz]/2, where z is the 
atomic position along the axis Oz of light propagation 
and Uo is the depth of the optical potential. The damping 
arises from stochastic transitions between the two optical 
potentials U±. These stochastic transitions correspond 
to the absorption and subsequent spontaneous emission 
of photons. Quantitatively, the departure rates 7±^t( z ) 
from the |±) states can be written in terms of the photon 
scattering rate V as 7±^^ = T'(l ±cos2fcz)/9 |23. The 
level of dissipation can therefore be quantitatively char- 
acterized by the photon scattering rate T', which can be 
controlled experimentally by varying the lattice fields pa- 
rameters. 




7C 2% 

FIG. 1: Results of semiclassical Monte Carlo simulations 
for a sample of n = 10 4 atoms in an ac driven ID lin_Llin 
optical lattice. The average atomic velocity, in units of the 
atomic recoil velocity v r — hk/m, is shown as a function of the 
phase difference <f) between the two harmonics of the driving 
force, see Eq. (0. Different data sets correspond to different 
scattering rates V , expressed in units of the recoil angular 
frequency u) T . The lines are the best fits of the data with the 
function v/v r — Asm((f) — (f>o). The calculations were done 
for a lattice with depth Uo = W0E r . The parameters of the 
driving are: A — 1.5, B — 6, ao = 0.2 ■ n, uj = 0.75u;„, where 
uj v is the vibrational frequency of the atoms at the bottom of 
the wells. 
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We studied the atomic dynamics in the presence of 
bi-harmonic driving by semi-classical Monte-Carlo simu- 
lations |22 |. For a given optical potential depth Uq, we 
calculated the average atomic velocity v as a function of 
the phase difference 4> between driving fields, for different 
values of the scattering rate. The results of our calcula- 
tion are shown in Fig. ^ The atomic current shows a 
dependence of the type of Eq. [21 with the phase lag (f>o 
vanishing in the Hamiltonian limit (r' — > 0) and increas- 
ing for increasing scattering rate. Values for 4>$ as a func- 
tion of the scattering rate, determined by fitting data as 
those of Fig.^with Asin(</> — <fio), are reported in Fig. [21 
for two different values of the driving force amplitude. 
It can be seen that although the magnitude and sign of 
the phase lag </>o are different for the two different driv- 
ing strengths considered, the general behavior described 
above is observed in both cases. For completeness, we 
mention here that also the sign and magnitude of the 
amplitude A varies depending on the driving strength. 
We notice that the dependence of the sign of 4>o and A 
on the driving strength is consistant with previous obser- 
vations of current reversals as a function of the driving 
amplitude at a g iven relative phase cj) between the two 
driving fields (24J . 

These results for the phase lag (/>n_are in agreement with 
general symmetry considerations .14] and also extend the 
validity of the theoretical results obtained for an ensem- 
ble of interacting particles 0] to the present system of 
non-interacting atoms, with the dissipation associated to 
the scattering of photons. 
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FIG. 2: Numerically calculated phase shift 4>o as a function of 
the scattering rate V . The phase shift is determined by fitting 
data as those in Fig.^with the function v/v r — Asin(<f) — cj>o). 
Different data sets correspond to different driving strength, 
i.e. to different values of A, with the ratio of the strengths 
of the two harmonics kept fixed: B — 4 A. All the other 
parameters are the same as for Fig. Q 

The experimental procedure followed closely the ap- 
proach of our numerical simulations. Pre-cooled cesium 
atoms are loaded in a near-resonant optical lattice. The 



phase modulation a(t) is then slowly turned on. By di- 
rect imaging the atomic cloud with a charge-coupled de- 
vice camera, we then derived the average atomic veloc- 
ity. By repeating the experiment for different values of 
the phase difference 0, we determined the average atomic 
velocity as a function of the phase <f>. 

Different sets of measurements were taken for different 
values of the scattering rate V at a constant depth of 
the optical potential. This was done by varying simul- 
tanously the intensity II and detuning A of the lattice 
beams, so to keep the potential depth Uq cx /l/A con- 
stant while varying the scattering rate V cx Il/A 2 - We 
notice that as II and A can be varied only within a fi- 
nite range, we cannot suppress completely dissipation, 
i.e. obtain V = 0. However, as we will see, for the driv- 
ing strength considered in the experiment, the smallest 
accessible scattering rate results in a phase shift which is 
zero within the experimental error, i.e. this choice of pa- 
rameters well approximates the dissipationless case. By 
then increasing V it is possible to investigate the effects 
of dissipation. 
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FIG. 3: Experimental results for the average atomic veloc- 
ity, in units of the recoil velocity, as a function of the phase 
4>. The lines are the best fit of the data with the function 
v/v r — Asin((j> — 4>o). The optical potential is the same for 
all measurements, and corresponds to a vibrational frequency 
uj v = 2n ■ 170 kHz. Different data sets correspond to different 
lattice detuning A, i.e. to different scattering rates as the 
optical potential is kept constant. The data are labeled by 
the quantity F s = [uj v / (2iv)] 2 / A which is proportional to the 
scattering rate. The parameters of the driving are uj = 2ir- 100 
kHz, A = 1, B = 4, Qo = 27.2rad. The errors on v/v r are of 
the order of 0.05. 

The results of our measurements, reported in Fig. |3 
demonstrate clearly the phenomenon of dissipation- 
induced symmetry-breaking. In agreement with our nu- 
merical calculations and with previous theoretical work 
[Til 0] , the measured current of atoms is well approx- 
imated by Asm((j) — cf>o). Therefore, by fitting data 
as those reported in Fig. |3| with the function v/v r = 
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Asm(4>— 4>q) we were able to determine the phase shift 4>q 
as a function of T', as reported in Fig. 0] The measured 
phase shift 0o is zero, within the experimental error, for 
the smallest scattering rate examined in the experiment. 
In this case, no current is generated for <fi — nn, with n 
integer, as for this value of the phase the system is invari- 
ant under time-reversal transformation. The magnitude 
of the phase shift (/>o increases at increasing scattering 
rate, and differs significantly from zero. The nonzero 
phase shift corresponds to current generation for (f> = mr, 
i.e. when the system Hamiltonian is invariant under the 
time-reversal transformation. This clearly demonstrates 
the breaking of the system symmetry by dissipation. 
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FIG. 4: Experimental results for the phase shift <j)o as a func- 
tion of F s = [lu v I (2ty)] 2 I A, which is proportional to the scat- 
tering rate. All the other parameters are kept constant, and 
are the same as for Fig. 



In conclusion, we demonstrated experimentally the 
phenomenon of dissipation-induced symmetry breaking 
with cold atoms in an optical lattice. We analyzed 
the atomic dynamics in an ac driven periodic optical 
potential which is symmetric in both time and space. 
These symmetries forbid the generation of a current. We 
showed that in the presence of dissipation the symmetry 
is broken, and a current of atoms through the optical lat- 
tice is generated as a result. Our results also show the 
generality of the phenomenon, in particular extending the 
results |l5| previously obtained for an ensemble of inter- 
acting particles in the specific framework of the kinetic 
Boltzmann equation to a system in which the dissipa- 
tion mechanism is of completely different nature. The 
present work is of relevance for the research on control of 
transport by time-dependent fields in a variety of system, 
ranging from optical tweezer set-ups [27j to quantum dots 
and wires 281. 



We thank EPSRC, UK and the Royal Society for fi- 
nancial support. 



[1] A. Ajdari and J. Prost, C.R. Acad. Sci. Paris 315, 1635 
(1992). 

[2] M. Magnasco, Phys. Rev. Lett. 71, 1477 (1993); R. Bar- 
tussek, P. Hanggi, and J. P. Kissner, Europhys. Lett. 28, 
459 (1994). 

[3] A. Ajdari, D. Mukamel, L. Peliti, and J. Prost, J. Phys. 

I (Paris) 4, 1551 (1994). 
[4] M.M. Millonas and M.I. Dykman, Phys. Lett. A 185, 65 

(1994). 

[5] M.C. Mahato and A.M. Jayannavar, Phys. Lett. A 209, 
21 (1995). 

[6] DR. Chialvo and M.M. Millonnas, Phys. Lett. A 209, 
26 (1995). 

[7] For a review on Brownian motors see R.D. Astumian, 
Science 276, 917 (1997); R.D. Astumian and P. Hanggi, 
Phys. Today 55, 33 (2002); P. Reimann, Phys. Rep. 361, 
57 (2002); P. Reimann and P. Hanggi, Appl. Phys. A 75, 
169 (2002); see also P.V.E. McClintock, Nature 401, 23 
(1999). 

[8] F. Jiilicher A. Ajdari, and J. Prost, Rev. Mod. Phys. 69, 
1269 (1997). 

[9] F. Marchesoni, Phys. Rev. Lett. 77, 2364 (1996). 
[10] M.I. Dykman, H. Rabitz, V.N. Smelyanskiy, and 

B.E. Vugmeister, Phys. Rev. Lett. 79, 1178 (1997). 
[11] I. Goychuk and P. Hanggi, Europhys. Lett. 43, 503 
(1998). 

[12] Ya.M. Blanter and M. Biittiker, Phys. Rev. Lett. 81, 



4040 (1998). 

[13] D.G. Luchinsky, M.J. Greenall and P.V.E. McClintock, 

Phys. Lett. A 273, 316 (2000). 
[14] S. Flach, O. Yevtushenko, and Y. Zolotaryuk, Phys. Rev. 

Lett. 84, 2358 (2000). 
[15] O. Yevtushenko, S. Flach, Y. Zolotaryuk and 

A. A. Ovchinnikov, Europhys. Lett. 54, 141 (2001). 
[16] P. Reimann, Phys. Rev. Lett. 86, 4992 (2001). 
[17] C. van der Broeck, R. Kawai, and P. Meurs, Phys. Rev. 

Lett. 93, 090601 (2004). 
[18] H. Linke, W. Sheng, A. Lofgren, H. Xu, P. Omling, and 

PE. Lindelof, Europhys. Lett. 44, 341 (1998). 
[19] H. Linke, T.E. Humphrey, Lofgren, A.O. Sushkov, 

R. Newbury, R.P. Taylor, P. Omling, Science 286, 2314 

(1999). 

[20] S. Weiss, D. Koelle, J. Miiller, R. Gross, and K. Barthek, 

Europhys. Lett. 51, 499 (2000). 
[21] J.E. Villegas, S. Save'ev, F. Nori, E.M. Gonzalez, 

J.V. Anguita, R. Garcia, J.L. Vincent, Science 302, 1188 

(2003). 

[22] For a recent review of optical lattices, see G. Grynberg 
and C. Mennerat-Robilliard, Phys. Rep. 355, 335 (2001). 

[23] M. Schiavoni, L. Sanchez-Palencia, F. Renzoni and 
G. Grynberg, Phys. Rev. Lett. 90, 094101 (2003). 

[24] P.H. Jones, M. Goonasekera, and F. Renzoni, Phys. Rev. 
Lett. 93, 073904 (2004). 

[25] R. Gommers, P. Douglas, S. Bergamini, M. Goonasekera, 



■5 



P.H. Jones and F. Renzoni, Phys. Rev. Lett. 94, 143001 
(2005). 

[26] F.-R. Carminati, M. Schiavoni, L. Sanchez-Palencia, 
F. Renzoni and G. Grynberg, Eur. Phys. J. D 17, 249 
(2001). 



[27] S.-H. Lee, K. Ladavac, M. Polin, and D.G. Grier, Phys. 

Rev. Lett. 94, 110601 (2005). 
[28] S. Kohler, J. Lehmann, P. Hanggi, Phys. Rep. 406, 379 

(2005). 



